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ABSTRACT 

Interacting galaxies often have complexes of hundreds of young stellar clusters of individual masses 
~ 10 4-6 M in regions that are a few hundred parsecs across. These cluster complexes interact dynamically, 
and their coalescence is a candidate for the origin of some ultracompact dwarf galaxies (UCDs). Individual 
clusters with short relaxation times are candidates for the production of intermediate-mass black holes of a few 
hundred solar masses, via runaway stellar collisions prior to the first supernovae in a cluster. It is therefore 
possible that a cluster complex hosts multiple intermediate-mass black holes that may be ejected from their 
individual clusters due to mergers or binary processes, but bound to the complex as a whole. Here we explore 
the dynamical interaction between initially free-flying massive black holes and clusters in an evolving cluster 
complex. We find that, after hitting some clusters, it is plausible that the massive black hole will be captured 
in an ultracompact dwarf forming near the center of the complex. In the process, the hole typically triggers 
electromagnetic flares via stellar disruptions, and is also likely to be a prominent source of gravitational radi- 
ation for the advanced ground-based detectors LIGO and VIRGO. We also discuss other implications of this 
scenario, notably that the central black hole could be considerably larger than expected in other formation 
scenarios for ultracompact dwarfs. 



1. INTRODUCTION 

Several bound systems of young, massive clusters in col- 
liding galaxies have been observed using the Hubble Space 
Telescope (HST). The best studied case is the Antennas galax- 
ies (NGC 4038/4039), the nearest examp le of two colliding 
disc galaxies listed in the Toomre ] ( 1977 ) sequence. HST 
observations reveal in this system the existence of relatively 
small regions (compared with the size of the galaxies) har- 
bouring hundreds or thousands of y oung clusters ( Whi tmore l 
19991 |Whitmore|[2006l ). In particular, Whit more| 



et al. 



2010 



et al.| ( 2010 ) observed 18 areas ("knots") of sizes spanning 
100-600 pc which contain hundreds of clusters. The mass 
function of those systems, which we will henceforth refer to 
as "Cluster Complexes" (CCs), is 



dN/dM oc M^, 



(1) 



with f3 = -2. 10 ± 2 in the range M ~ 10 4 ~ 5 M© (see also 
Zhang & Fall|1999| ). |Bastian et aL] ( |2006l ) found in the same 
system low-mass CCs with masses around 1O 6 M and diam- 
eters of some 100-200 pc. One of the best studied CCs in the 
Antennae galaxy is "knot S", with a total mass of 1O 8 M and a 



total radius of ~ 450pc (Wh itmore et al . 1 1 999] ) . Othe r galaxies 
with recent ly disc overed CCs include NGC 7673 (|Homeier 



eTaL][2002l ), M82 ( [Konstantopoulos et al.|[2QQ9] ), NGC 6745 
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Ode G rijs et al.||200 3|), Stephan's Quint et ((Gallagher et al. 
2001 ) and NGC 922 ( |Pellerin et a l. 2010 ). 

CCs are bound systems (|Kroupa|199 8 ; Fellhauer & Kroupa 
[20051 IBruens et al.|20lT|| Whitmore et al.|2010| ) and on rela- 
tively short time-scales at least som e of their memb er cl usters 
will merge to form a si ngle object. Kroupa (1998 ) and |Fell- 
|hauer & K roupa (2005 ) have postulated CCs as the breeding 
ground of Ultra-Compact Dwarf Galaxies (UCDs). Following 
this idea, Brue ns et"aL] ( |201l| ) performed A/-body simulations 
of CCs with different total masses (10 5 5 - 10 8 M ) and initial 
Plummer radii 10- 160 pc. They conclude that UCDs, Ex- 
tended Clusters (ECs) or even large Globular Clusters (GCs) 
might be the product of an agglomeration of clusters in CCs. 
They find in their simulations that almost all members of a 
CC merge in less than 1 Gyr. In some cases this timescale can 
be as short as 10 Myr. By the end of their simulations a very 
massive cluster forms in the centre of the CC, with a mass of 
26-97% the mass of the initial CC and a radius of ~ 50 pc. 

Theoretical and numerical studies show that at least a frac- 
tion of young star clusters could host intermediate-mass black 
holes (IMBHs, black holes with masses ranging between 
10 2-4 M ) at their centres. A possible formation path is that in 
a young cluster, the most massive stars sink to the centre due 
to mass segregation. After a high-density stellar region forms, 
stars start to collide and merge with each other. A number of 
numerical studies with rather different approaches show that, 
under these circumstances, at least one of the stars increases 
its mass rapidly in a process of runaway collisions (Portegies 



Zwart & McMillan 2000; Gurkan et al. 2004; Portegies Zwart 



et al.||2004| IFreitag et al.||2006a|b ) . Nonetheless, there are 
a number of open questions regarding this process. One of 
the main uncertainties is the role of stellar winds. In prin- 
ciple at approximately solar metallicity winds may limit the 
mass of this very massive star (VMS) to a few hundreds o f 
solar masses rather than a few thousands (Belkus et al.|2007] ). 
Nevertheless we note that this requires a substantial extrapo- 
lation of already uncertain wind loss rates to stellar masses an 



2 



order of magnitude beyond what is observed. Also, the col- 
lision process might lead to lumpy bags of stellar cores in an 
extended envelope rather than to relaxed stars near the end of 
the runaway co llision (M. Davies , private communication). In 
addition, when |Suzuki et al.| ( |2007| ) combined direct A/-body 
simulations with smooth particle hydrodynamics (SPH) they 
found that stellar winds would not hinder the formation of the 
VMS. It is thus possible but not certain that IMBHs can form 
in young clusters. We will assume their existence as a work- 
ing hypothesis. 

Apart from the obvious interesting implications for mod- 
els of galaxy formation and, in particular, of UCDs, merg- 
ers of clusters in CCs are a powerful source of gravitational 
waves if these harbour central IMBHs in their respective cen- 



tres ( Amaro-Seoa ne & Fre itag 2006; Amaro -Seoane et al. 
|2010 , 2009; Amaro-Seoane & Santamarf a|2010| ). In particu- 
lar, Amaro-Seoane & Freitag ( 2006 ) showed that such merg- 
ers would lead to the formation of an IMBH binary, which 
would merge in a time scale as short as ~ 7 Myr. Such a 
merger would be easily detected with space-borne observato- 
ries and also with ground-based detectors such as Advanced 
LIGO or Advanced VIRGO (AdLIG O/AdVIRGO) if it oc- 
curs within ~ 2 Gpc (Frege au et al.||2006|). Using more re- 



alistic waveforms including spins, Amaro-Seoane & Santa 



mana 



( |2010| ) find that the detection distance is increased sig 
nificantly, up to an orientation- averaged distance of ~ 5 - 12 
Gpc, depending on the spin configuration and mass ratios. In 
the case of the Einstein Telescope (ET), the same authors find 
that the maximum redshifts for ET are z ~ 10, which implies 
that binaries of IMBHs will be a cosmological probe. 

Numerical relativity simulations show that during the 
merger of the holes, gravitational radiation is emitted asym- 
metrically with the size of asymmetry depending on the mass 
ratio of the tw o black holes and on their spin magnitude 
and orientation ([Zlochower et al.|2010[|Gonzalez et al.|2007 ; 



Lousto & Zlochower 2008; Bak er et al.||2008i |Lousto et al. l 
20101 ICampanelli et , al.|2007a|b[|Heary et al. 12009 UPretorius 
2005; Sopuert a"eTal.||2006[ |Herrmann et al.|l2007a|b[fB oyle 
& Kesden|2008Hvan Meter et al.|2010||Lousto "STZ lochower 



2011b ) If this recoiling velocity exceeds a few times the ve- 



locity dispersion of the merged cluster, then the IMBH leaves 
the host cluster. There is a massive black hole at large in 
the CC. Even if an IMBH escapes from one cluster, it might 
still be bound to the CC as a whole, which means that it has 
the possibility of interacting with other clusters and, perhaps, 
their IMBHs. 

In this article we address the formation of ultra-compact 
dwarf galaxies by the agglomeration of young clusters in CCs, 
along with the role of one or more recoiling IMBHs, using 
direct- summation A/-body simulations. For this, we run a set 
of ~ 200 individual experiments in which we vary mass ra- 
tios, relative speeds, and impact parameter to study in detail 
the interaction between a single IMBH and a cluster. We then 
study the interaction of one or more IMBHs at large in a CC 
with individual clusters with an additional set of A/-body sim- 
ulations. We correct for the trajectory of the IMBH, based on 
point dynamics and the mass loss in the individual clusters, 
by using the previous 200 experiments. We also follow the 
growth of a seed UCD in a CC and record all stellar disrup- 
tions triggered by the presence of the IMBH(s). For realistic 
models of CCs we find that the IMBH(s) end up captured by 
the seed UCD or by a smaller cluster which is close to the 
UCD. Thus, if the fraction of IMBHs in the CC (/. from now 
onwards) is not zero, this is a process of allocating one or 




FIG. 1. — Geometry for the initial conditions of the parabolic collision, 
in the COM of the IMBH-cluster system. To obtain the grid displayed in 
Fig.j2j, we systematically vary ^min? the relative velocity and the mass ratio 
between the IMBH and the cluster. 

more IMBH at the very centre of a UCD. 

2. INTERACTIONS BETWEEN A RECOILING IMBH AND AN 
INDIVIDUAL YOUNG CLUSTER 

In this section we make a study of the parameter space for a 
collision between a recoiling IMBH and an individual young 
cluster in a CC. We run a set of ~ 200 direct A/-body simula- 
tions to build a grid which we will later use in our simulations 
of the IMBH in the CC, as explained in the introduction. Ini- 
tially we set the IMBH and the cluster on an orbit with positive 
relative speed and thus positive total e nergy in the initial state , 



i.e. a hyperbolic orbit, as described in |Amaro-Seoane| ([2006 ). 
We schematically show this for reference in Fig. M and fol- 
low a similar notation. The initial trajectory of the IMBH 
would bring it within a minimum distance d m { n of the cluster 
centre if the cluster was replaced by a point particle. In the 
centre-of-mass reference frame (COM), 



x. = A d d, 
x c i = -A. d, 

V. = A c iV r el, 

v c i = -A. v rd (2) 

where v re i is the relative velocity of the two objects d is their 
separation vector, x. c i are the positions of their centres, and 
A.,ci = m, , c i/(M 9 +M c \). 

The number of stars in the cluster is always N^ = 3x 10 4 and 
we use for the ir initial distribution a King model of concen- 
tration Wo = 7 ( |King|1966j|Peterson & King|1975| ). The sim- 
ulations are performed with the direct- summation NBODY4 
programme of | Aarseth] ( [2003] ). This choice was made for the 
sake of the accuracy of the study of the orbital parameter evo- 
lution of the IMBH and mass loss in the cluster; this numer- 
ical tool include s both the KS regularisation ( Kustaanheimo 
|& Stief el]| 1 965| > and chain regularisation, which means that 
when two or more particles are tightly bound to each other or 
the separation is very small during a hyperbolic encounter, the 
system becomes a candidate to be regularised in order to avoid 
problematical small individual time steps. The basis of direct 
Af-body cod es relies on a Hermite integrator scheme ( Aarseth 



|1999| |2003| ) for which we need not only the accelerations but 
also their time derivatives. This extra computational overhead 
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is necessary for us to follow reliably the orbital evolution of 
every single star (or IMBH) in our system. 

At the end of an Af-body run, we need to identify the par- 
ticles that are still forming a bound cluster, the particles that 
are bound to the IMBH, and the particles that have become 
unbound. We also need to know whether the IMBH has been 
captured by the stellar cluster. We have therefore developed 
an iterative algorithm. To initialize the procedure, we make a 
(computationally) quick guess of which particles are bound to 
the cluster and which ones form a bound group including the 
IMBH (the "IMBH group"). Note that a given particle can be 
in both groups, for instance if the IMBH has been captured by 
the cluster and has sunk to its centre or is orbiting it. For this 
first guess, stellar particles are considered bound to the IMBH 
group if they are bound to the IMBH (i.e., we do not take into 
account the self-gravity of the bound stars themselves). 

For the first-guess cluster, one assumes that its centre cor- 
responds to the median position of all stellar particles, i.e. the 
x, y and z components of the "centre" are taken to be the me- 
dians of the corresponding components of the positions of all 
the stellar particles. The median turns out to be a much more 
robust estimate of where the bulk of the particles is, com- 
pared to the average position or the centre of mass (i.e. the 
mass-weighted average position) as those quantities are very 
sensitive to the the positions of a few particles ejected at large 
distances from the rest. For this first guess, the 90 % of the 
particles closest to this median position are assumed to be part 
of the cluster. 

For the first iteration, we have to compute the binding en- 
ergy of a particle relative to the cluster group, hence we need 
to know the velocity of that group. To estimate the velocity of 
the cluster in the first-guess attribution, we take the average 
velocity of the 10% of the particles closest to the assumed 
centre. This number is sufficient to avoid large fluctuations 
due to individual particle velocities ("random velocities"). On 
the other hand, taking a significantly larger fraction of parti- 
cles is neither necessary nor advisable as it is not yet known 
which particles are actually bound together as a cluster. We 
have to make sure that the velocity defined in the procedure 
is a good estimate of that of the actual bound cluster. Other- 
wise, the kinetic energies relative to this first-guess cluster are 
biased towards high values and the iterative procedure fails 
at identifying a bound cluster. The iterations proceed as fol- 
lows: For each particle, the binding energies relative to the 
cluster and the IMBH group are computed. For this, we have 
to estimate the position of the centre of each group and its ve- 
locity. For the IMBH group, they are fixed to the values of the 
IMBH itself. For the cluster, the centre position and velocity 
are defined to be the mass-weighted mean values for all parti- 
cles within half a "typical size" of the previous estimate of the 
centre. The typical size of the cluster is the harmonic mean of 
the distance to its centre (for all particles considered bound to 
it): 

^typ=^harm=M cl ^^^ . (3) 

One advantage of defining R typ using the harmonic mean, in- 
stead of using the half-mass radius or some other Lagrangian 
radius, is that this does not require a sorting of the particles. 
The gravitational energy is computed assuming a spherical 
mass distribution, i.e., as if each particle bound to a group 
(cluster or IMBH group) was a spherical shell of matter, of ra- 
dius Ri centred on either the IMBH position or the estimated 



centre of the cluster. Typically, the attributions of the par- 
ticles to either or both groups converge after fewer than ten 
iterations. 

At the end, the attributions are cleaned up in the following 
way. If a stellar particle belongs to both the cluster and the 
IMBH group, the binding energies to both structures are com- 
pared. It will be kept as member of the IMBH group only if 
the binding energy to the IMBH group is larger than to the 
cluster group. In that case, it will also be kept as member 
of the cluster only if the IMBH itself is bound to the cluster. 
This reduces the number of double-members in a reasonable 
way, still allowing for situations such as the IMBH having 
captured some stars while being itself on a bound orbit around 
the (main) cluster. 

Finally, to interpret the results, we allow for three different 
outcomes. A merger is when the IMBH group is bound to the 
cluster (as determined assuming each group is a point mass) 
and the distance between the centres of the groups is smaller 
than the sum of the ^ typ 's. A satellite situation arises when the 
two groups are bound but the distance between their centre is 
larger than twice the sum of the ^ typ 's. A fly by is when the 
groups are unbound and the distance between their centres is 
larger than either the sum of the total extent of each group or 
five times the sum of the ^ typ 's. Any other situation would 
be considered as unknown but does not occur if the A/-body 
simulation has been carried out for a sufficient duration. 
In Fig.(|3]) and [4] we show two particular cases for the IMBH - 
cluster interaction in the COM frame which, although not rep- 
resentative for the whole sample displayed in Fig. ([2]), are in- 
teresting in terms of the dynamics of the system|j. In the first 
case J m i n = 1 , which leads to an almost head-on collision be- 
tween the IMBH and the cluster. Nonetheless, because of the 
low relative velocity and mass ratio, the interaction does not 
lead to a huge mass loss from the cluster. Even if at T = 45.60 
Myr the IMBH and cluster seen to be unbound, the IMBH is 
still forming a binary with the COM of the cluster and, hence, 
the semi-major axis decays again. After some 154 Myrs the 
IMBH settles down to the centre and is captured. In the sec- 
ond figure, the larger mass ratio has a significant impact in 
terms of mass loss. Already after 11.62 Myr the IMBH has 
captured some stars from the cluster, which remain bound to 
the trajectory of the hole and follow its trajectory. This satel- 
lite and the IMBH are nevertheless still gravitationally bound 
to the cluster and hence fall back again. The higher mass in 
the IMBH-satellite system leads to a rather large mass loss 
from the original cluster. After 80.50 Myrs the IMBH is at 
the centre of the remaining cluster. 

3. INTERACTIONS BETWEEN A SINGLE RECOILING IMBH AND 
CLUSTERS IN A CC 

3.1. Integrator and first considerations 

Now that we have completed the grid of individual IMBH- 
cluster interactions, we can explore the scenario in which one 
IMBH is at large in a CC, interacting with many different 
IMBH on its way, either to an eventual escape from the CC, 

8 The interested reader can visit 

I http : / /members . aei .'m pg . de/amaro- seoane/ | 
|ul tr a- compact -dwarf- gal axies| 

for movies based on the results of the figures (the last URL is a 3D version 
of the second movie). The encoding of the movies is the free OGG Theora 
format and should stream automatically with a gecko-based browser (such as 
mozilla or firefox) or with chromium or opera. Otherwise please see e.g. 

http : / /en . wikipedia . org/wiki/Wikipedia : Media_help_ 

l(Qgg)l 

for an explanation on how to play the movies. 
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Fig. 2. — Outcomes of all 196 simulations of encounters between a cluster with King parameter Wo = 7 and an IMBH. Each panel shows the results for a given 
mass ratio A4 m /A4 c \. The abscissa of each plot is the minimum distance d^, computed assuming 2-body dynamics, in units of the half-mass radius R^. The 
ordinate is the relative velocity at infinity Voo, in units of Vh = (GM. C \/Rh) 1 ^ 2 , a typical velocity dispersion for the cluster. Solid round dots show "mergers", i.e., 
cases where the IMBH has been captured by the cluster and has settled at its centre. Solid triangles are cases in which the IMBH is orbiting the cluster (a merger 
is likely to be the long-term outcome). Open squares are "fly-throughs". The number just below a symbol (in blue in the on-line colour version) is fractional mass 
loss from the cluster in percent. The second, lower number (in orange in the on-line colour version) is the fractional reduction in specific binding energy of the 
cluster, also in percent. A number above a symbol indicates how many stellar particles are bound to the IMBH (when it has not merged with the cluster). 

or down to the very centre, where the seed of a UCD is form- 
ing by the mergers of clusters. In this section we will assume 
/. = 2/Ntot, where N tot is the total number of clusters in the 
CC; that is, only two clusters in the whole CC harbour an 
IMBH and we also assume that they have coalesced and the 
merged hole escaped from the host cluster. As we will see, 
the presence of the IMBH triggers stellar disruptions in indi- 
vidual clusters of the CC, which could potentially represent 
a fingerprint of this process. In next section we will look at 
larger values of /. . 
The numerical code that we us e for the simulations of the 

Kokko-| 



CC and t he IMBH is Myriad ( |Konstantinidis & 
tas 2010), which uses the Hermite fourth-order predictor 



corrector scheme with block time steps (Makino & Aarseth 
|1992| ) for advancing the particles in time, while the accel 
erations and their d erivatives are computed using GRAPE-6 
(Mak ino et al.||2003] ) special purpose computers. Close en- 
counters between stars are detected using the GRAPE-6 and 
evolved with a time- symmetric Hermite fourth-order integra- 



tor ( [Kokubo et al.|1998) . Even though the code was originally 
designed for dynamical simulations of stars in star clusters, 
its flexible modularity made it easy to adapt to our particular 
problem. In particular, we assigned a radius to each particle 
representing a cluster, and we allowed clusters to merge with 
each other whenever the distance was smaller than the sum of 
the radii. In the simulations the IMBH is also a particle with 
a radius set to its Schwarzschild radius. 

From the individual IMBH-cluster simulations presented 
previously we have data for the outcomes based on the mass 
ratio Mm/Mch the distance of closest approach between 
IMBH and cluster, and the relative velocity of the two ob- 
jects and, thus, the change in kinetic energy of the IMBH. We 
use these results to correct the position and velocity of the 
IMBH after each interaction with a cluster in the simulation 
of the CC. This also provides us with information about the 
number of stellar disruptions triggered by the IMBH, as well 
as the characteristics of the cluster which captures the IMBH 
(if any). If a capture does occur, the simulation finishes and 
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Fig. 3. — Projection in the X-Y plane of all trajectories of the stars (star symbols) in a cluster and the IMBH (red circle) for 12 different moments in the 
interaction. In this particular case, the process leads to the capture of the IMBH. For visibility, the radius of the IMBH and the stars has been artificially magnified. 
We also depict the previous 60 positions of the IMBH with a solid, green line. The mass ratio between the IMBH and the cluster is 0.01, the minimum distance 
of approach of the COM of the cluster and the IMBH is fimin — 1 and Vqo — 1 km s 

then we record the position of the "trapping" cluster in the CC. 
Another possible termination of the simulation is if the IMBH 
leaves the CC, because its speed is high enough to escape the 
complex. 

3.2. Assumptions for the initial conditions of the CC and the 

IMBH 

Initially we fix the radius of the CC, 7?cc> to a typical value 
coming from observational data and populate it with individ- 



ual clusters following equation [T] In particular, in the "knots" 
of the Antennas galaxy one observes a mass distribution with 
n = -2. The number of observed individual clusters in CCs is 
of the order of 100, but the actual number might actually be 
thousands, most of which are simply too faint t o be observed 
(as discussed in e.g. Fellhauer & Kroupa 2005 ). We set the 
total mass of the CC to a typical observed value, M C c ~ 10 6 - 
1O 8 M . The individual clusters have half-mass radii ranging 
between 0.5 and 4 pc and are distributed initially in the CC 
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Fig. 4. — Same as in Fig.jiJ for 15 different times. The mass ratio in this case is 0.333, d m \ n = 5 and Vqo = 3kms 1 . 
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Fig. 5. — Kinetic energy difference between the initial and final kinetic 
energy (7]dnf, i and f respectively) normalized to the initial energy for all 
collisions between the IMBH and the clusters resulting in a fly-through for 
all TV-body simulations. 



following a Plummer model (Plumme r|1911| ) with a cut off ra- 
dius (see table [T). The masses of the clusters are discrete and 
come from the M./ M c \ ratios that were used in the IMBH- 
cluster A/-body simulations. Then, for M. = 5 x 1O 3 M 
and M. / M c \ = 0.01 , 0.033 , 0.1,0.33,1,2, the discrete masses 
of the clusters in the CC are 5 x 10 5 M o , 1.51515 x 
1O 5 M ,5 x 1O 4 M , 1.51515 x 1O 4 M ,5 x 1O 3 M and2.5 x 
1O 3 M . When two clusters collide in the CC simulation, we 
assume a 20% mass loss, so the cluster product of the merger 
of two individual clusters has a mass which is 80% of the sum 
of the masses and a new radius, equal to the radius of the more 
massive cluster plus the 20% of the sum of the radii of the two 
clusters. 

The IMBH in the merged cluster is assumed to be the prod- 
uct of a merger of two IMBHs that were located at the centres 
of two merging star clusters. We assume that this happened 
close to the centre of the CC, where most of the individual 
cluster-cluster collisions take place, because this is where the 
numerical density of clusters is highest. Hence, we initially 
place the IMBH at the centre of the CC. We choose a mass of 
A4 9 = 5 x 10 3 Mq, which determines the masses of individ- 
ual clusters from the grid given in the previous section. The 
recoil speed of t he merged IMBH could in principle be up to 
5000 km s" 1 ((Herrmann et al.||2007a|b[ Boyle & Kesden 
TI20TTBT 



2008[ |Lousto & Zlochower 201 l b|a| )for optimal mass ratios, 
spins, and spin orientations. The recoils of greatest interest 
to our present study are in the ~ 100 km s -1 range, because 
the merged IMBH will then escape from its host cluster but 
be bound to the CC as a whole. It is difficult to judge how 
representative this will be for the mergers of actual IMBHs 
in CCs. Assuming the spin orientations are random, speeds 
in this range are characteristic of mass ratios q ~ 0. 1 for sub- 
stantial sp ins, or spins a/M ~ 0.1 for mass ratios co mparable 



to unity ( |Lousto et al.|2010[|van Meter et al.|2010| ). For our 
purposes we will study the case of v reco ii = 100 km s -1 . At this 
speed, the escape time from a cluster of total radius ~ lOpc is 
~ O.lMyr. Hence we simply place the IMBH initially at the 
centre of the CC, not bound to any cluster, and assume that it 
recoils in a random direction. 

For the evolution of the recoiling IMBH we must take into 
account the loss of kinetic energy every time it hits a cluster. 
In figure ([5]) we can see the distribution of the resulting kinetic 



energy after a hit for all fly-by simulations of figure ([2]). We 
can clearly see that a bit more than 50% of the systemsnave a 
loss of about ~ 5% after a hit. This loss of energy will result 
into a rather negligible deacceleration of the IMBH, so that 
it will have more chances to escape the CC, and it will also 
lead to a lower number of stellar tidal disruptions. On the 
other hand, a bit less than 50% of all "fly-throughs" have at 
least over ~ 5% of relative loss after one hit. This situation 
is more appealing from a dynamical standpoint, and therefore 
we will first address it. In the next sections we will assume an 
average loss of 20% for the "fly-throughs" hits, and in section 
[5] we will briefly explore the other situation. 

Our parameter space consists of the number of clusters N 
and the initial radius of the CC, Rqc- The total mass M C c 
of the CC is a consequence of N, because the masses of the 
clusters are assumed to follow a power law. The total radius 
that we use varies from 45 pc to 330 pc. Given the mass and 
the size of the CCs, the initial escape speeds at the centres 
of the CCs are between 27- 137kms _1 . All details for all 
simulations are given in table [T] 

3.3. Results of the simulations 

In figure [8] we present the results of our 34 CC simulations. 
In simulations A1-A5, B1-B5, C1-C4, D1-D4, E3-E4, F4, G4 
and H4 the IMBH escapes the CC after between zero and 
a few interactions with clusters. These cases correspond to 
smaller-mass CCs or to low initial concentrations. In simu- 
lations E1-E2, F1-F3, G1-G3, and H1-H4, which are more 
representative of observed CCs, the IMBH is captured in the 
CC after a significant number of interactions and ends up be- 
ing trapped by an individual cluster, which can be the UCD 
seed (cases El, E2, F2, F3, and H2). We show two particular 
cases which led to the capture of the IMBH in figures [6]and[7] 

The IMBH goes through a very large number of interactions 
with individual clusters until it is eventually trapped. This 
number depends on the density of clusters in the CC. In 6 
simulations, the IMBH gets captured by a cluster that has not 
yet merged with other clusters. In 5 simulations, the cluster 
that captures the IMBH is the central cluster of the CC, the 
seed UCD. We show in table [2] the details about the cluster 
that captures the IMBH, the distance from the centre where 
this takes place, and the mass of the most massive cluster in 
the system at the time of the IMBH-capture, i.e., the mass of 
the UCD seed. An interesting process in the dynamical evolu- 
tion of the system is that the IMBH triggers stellar collisions, 
i.e., stars are set on such an orbit that they collide and dis- 
appear from the system. We note that only in one case, in 
simulation Fl, one star was torn apart by the tidal forces of 
the IMBH acting on a star. The middle number next to each 
circle of figure [8] corresponds to star-star collisions triggered 
by the IMBH in the clusters. We can conclude that one should 
expect a star-star collision in a CC every 5 - 8Myr. In figure 
[9] we show the accumulated number of stellar collisions that 
led to a disruption in function of the time for simulation G3, 
as well as the accumulated number of hits between the IMBH 
and a cluster. 

The third number next to each circle of Figure[8]is the initial 
escape velocity at the centre of the CC. As it is obvious, CCs 
with values < 100km s -1 retain the IMBH due to our choice 
of the initial recoiling speed. An interesting case is simulation 
H4 in which the escape velocity is 84km s -1 , but the IMBH 
escapes because the system is initially not very concentrated 
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ID 


N 


Mcc(Mq) 


Rcc (pc) 




ID 


N 


M cc (M Q ) 


Rcc (pc) 


Al 


5x 10 2 


1.522 x 10 7 


45 




El 


3x 10 3 


4.32 x 10 7 


122 


A2 


5x 10 2 


1.522 x 10 7 


90 




E2 


4x 10 3 


5.75 x 10 7 


165 


A3 


5x 10 2 


1.522 x 10 7 


132 




E3 


4x 10 3 


5.75 x 10 7 


246 


A4 


5x 10 2 


1.522 x 10 7 


168 




E4 


4x 10 3 


5.75 x 10 7 


329 


A5 


5 x 10 2 


1.522 x 10 7 


255 












Bl 


1 x 10 3 


1.522 x 10 7 


90 




Fl 


5x 10 3 


7.18 x 10 7 


122 


B2 


1 x 10 3 


1.522 x 10 7 


128 




F2 


5x 10 3 


7.18 x 10 7 


165 


B3 


1 x 10 3 


1.522 x 10 7 


169 




F3 


5x 10 3 


7.18 x 10 7 


248 


B4 


1 x 10 3 


1.522 x 10 7 


252 




F4 


5x 10 3 


7.18 x 10 7 


330 


B5 


1 x 10 3 


1.522 x 10 7 


333 












CI 


2x 10 3 


2.9 x 10 7 


126 




Gl 


6x 10 3 


8.6 x 10 7 


122 


C2 


2x 10 3 


2.9 x 10 7 


167 




G2 


6x 10 3 


8.6 x 10 7 


165 


C3 


2 x 10 3 


2.9 x 10 7 


252 




G3 


6 x 10 3 


8.6 x 10 7 


248 


C4 


2x 10 3 


2.9 x 10 7 


336 




G4 


6x 10 3 


8.6 x 10 7 


330 


Dl 


3x 10 3 


4.32 x 10 7 


124 




HI 


8x 10 3 


1.14 x 10 8 


122 


D2 


3x 10 3 


4.32 x 10 7 


166 




H2 


8x 10 3 


1.14 x 10 8 


165 


D3 


3x 10 3 


4.32 x 10 7 


249 




H3 


8x 10 3 


1.14 x 10 8 


248 


D4 


3x 10 3 


4.32 x 10 7 


332 




H4 


8 x 10 3 


1.14 x 10 8 


330 



TABLE 1 

Simulation ID, number of clusters, total mass and cut-off radius of the CC. Note that the table is vertically split in two 

SUBTABLES. 



and the IMBH has only 2 interactions with clusters. In this 
case, the energy of the IMBH did not decrease enough to be 
trapped in the CC. Simulation G3 corresponds to the opposite 
situation. Even though the escape speed is the same as in H4, 
the IMBH remains in the system because the CC is denser, so 
that the IMBH has a chance of interacting significantly with 
clusters and, hence, of decreasing its kinetic energy below the 
threshold. In figure [10| we have the evolution of the velocity 
of the IMBH in simulation G3 compared with the escape ve- 
locity at the radius of the CC where the IMBH is. Initially, the 
escape velocity is lower than the velocity of the IMBH, ensur- 
ing the escape of the IMBH from the system, but the IMBH 
loses energy rapidly during the first few Myr because of its 
interactions with clusters. 

4. INTERACTIONS BETWEEN MULTIPLE RECOILING IMBHS AND 
CLUSTERS IN A CC 

In this section we investigate a scenario in which f. > I. 
We use the initial configuration of F3 as described in table [T] 
as our fiducial CC system and study the evolution of systems 
of five and ten IMBHs at large. For this, we set them initially 
close to the centre and allow them to be kicked off the host 
cluster at the same time, T = 0, as a simplifying assumption. 
In real systems there will be a time lag: 



7"bin — 7"run 7"IMBH 7"merg • 



(4) 



where r mn is the timescale for a cluster to evolve to the run- 
away phase, timbh is the timescale for the VMS to become 
unstable and form an IMBH and r merg is the timescale for 
the cluster to merge with another cluster. The phenomena in- 
volved are various and the assumptions inherent to r mn and 
7" imbh prevent realistic estim ates, as we explained in the intro- 
duction. On the other hand, Amaro-Seoane & Freitag ( 2006 ) 
estimate that T merg ~ 7 Myr, which compared to the timescale 
for the CC to reach the seed UCD phase, of the order of ~ 100 
Myr, is a rather short interval of time and can be regarded as 
instantaneous. In view of these arguments, we assume that the 
IMBHs are expelled instantaneously from their host clusters 
at different places of the CC at T = 0. 



In table [3] we show the results for the first simulation, in 
which we place five IMBHs around the centre, as indicated in 
column number three. IMBHs #2-5 have been distributed 
over the surface of a sphere of radius 17.32 pc and only one, 
#1, is very close to the centre, to avoid the artificial forma- 
tion of various binaries of IMBHs when we start the sim- 
ulation. We assign the holes initial recoil speeds between 
50- lOOkrns -1 and different directions and then let the sys- 
tem evolve. We find that after some ~ 34 Myr all IMBHs 
have been either captured by an individual cluster which is 
sinking the the centre due to DF, or formed a satellite with a 
cluster. In figure [TT] we show the CC at T = 62.37 Myr. We 
stop the simulation at that time because the satellites are con- 
suming all of the computational power. In the process and up 
to that time, there are 7 stars that have been disrupted in the 
CC, as we can see in the table. 

In table|4]we repeat the same exercise but for a system with 
10 IMBHs. The initial setup is identical to the previous one. 
We find that in this case three holes leave the system due to an 
increase in their kinetic energy. The rest of them have formed 
a hard binary with a cluster and will eventually be captured. 

5. LOWER KINETIC ENERGY LOSS 

In the simulations of the previous sections we assumed a 
loss of relative kinetic energy of ~ 20% for the hits that led 
to a fly-through, although it could be much larger than that, as 
we saw in figure [5] While this is true for a bit less than 50% 
of all systems, the rest of them had a peak in the distribution 
around ~ 5%. We have addressed the situation of a larger loss 
first, because it leads to more interesting effects from a pure 
dynamical standpoint. 

However, we deem it necessary to we repeat some exper- 
iments in the evolution of the CC to understand the other 
regime. Therefore, we repeat experiments Gl, G2, G3, G4, 
HI, H3, H4 of table [T] but this time we assume a loss of 5% 
after every hit for the fly-throughs. In figure ( [T2| ) we can see 
the results. We have reduced the exploration to the range of 
radii and total mass that could be more interesting for our 
analysis. We can see that although the total number of stel- 
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Fig. 6. — Formation of the UCD seed at the centre of the CC. We show a projection in the X-Y plane of all individual clusters for the simulation F3. The radii 
of the clusters have been artificially magnified, heavier members have larger sizes and darker colours relative to every panel for the sake of visibility. This means 
that even if the colours of the heaviest clusters in the last panel are as dark as the most massive ones in the first panel, the clusters in the last panel are heavier 
and larger. After 7.44 Mys we can already see how the more massive clusters start to agglomerate at the centre of the CC. Later, at T ~ 40 Myr, all of them are 
confined to the central part of the CC and in the last panel we can see that only a handful of clusters are heavy and a very massive cluster is sitting at the very 
centre, while lighter clusters occupy all of the remaining space. The mass of this very massive cluster is 2.9 x 1O 6 M0 and constitutes the seed of the UCD. See 
|http : / /members . aei .mpg . de/amaro- seoane /ultra- compact -dwarf- galaxies] model F3 for an animation of the process. 



lar disruptions is signficantly reduced, it is not zero. Also, in 
four configurations the IMBH at large is captured eventually 
by the forming CC. 

6. SUMMARY AND CONCLUSIONS 

In this work we have presented results that address the for- 
mation of UCD from young clusters, and the role of recoil- 
ing IMBHs in a CC. For that, we first ran a set of ~ 200 
direct- summation A/-body simulations that covers the param- 
eter space for individual IMBH-cluster encounters. We me- 
thodically varied the mass ratio between the IMBH and the 
cluster, the relative velocity, and the impact parameter. This 
allowed us to build a grid with the expected outcome of the 
interaction and the modification of the kinetic energy of the 
IMBH. Later we ran additional direct- summation A/-body 
simulations for a scenario in which one IMBH is at large in a 
CC. The IMBH is assumed to be the result of the coalescence 
of two holes, which led to the expulsion of the hole from the 
initial host cluster. We studied the dynamical evolution of 
this single IMBH in an evolving CC. Parallel to the individ- 



ual interactions between the IMBH and clusters in the CC, 
which are corrected using the above-mentioned table, clus- 
ters are colliding and merging with each other, which results 
in the formation of a run- away individual cluster, which typi- 
cally after ~ 100 Myrs contains almost all of the mass of the 
CC. This is what we designate "the seed of an ultra-compact 
dwarf galaxy", since this very massive cluster is the result of 
the successive amalgamation of smaller clusters in the initial 
distribution of the CC. 

We find that for realistic CCs (i.e. those which resemble 
observations, such as the knots of the Antennas), the IMBH 
is either eventually captured by the seed UCD (in those sim- 
ulations less dense initially) or by a smaller cluster (in the 
simulations with the largest concentrations of clusters at the 
centre) which, however, is close to the centre of the CC, so 
that it will in the course of time sink down to the very centre, 
where the seed UCD is settled. The typical timescale for this 
trapping is of about ~ 200 Myr. 

We can see this by estimating the dynamical friction time 
r DF . This is the timescale for the IMBH captured in a cluster 
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Fig. 7. — Same as in Fig. {6} but for simulation G3. In this case we show a zoom of diameter 600 pc. As in the first figure, after some ~ 100 Myr we have a 
very massive cluster at the centre and all other clusters are much lighter. The heaviest cluster at this time has a mass of 5.5 x 10 5 , while clusters with masses 
5.2 x 10 5 M^, 5.0 x 10 5 M^, 1.9 x 10 5 M^, 1.4 x 10 5 M^ an d 6.5 x 1O 4 M lie very close to the centre of the CC. Seej |http : / /members . aei . mpg . de/| 
|amar o- seoane/ ultra- compact -dwarf- gal axiesj model G3 for a movie of the figure. 



ID 


Coll 


T[Myr] 


Rcapttpc] 


M cl [M Q ] 


-MucdDVI©] 


r DF [Myr] 


TpntMyr] 


7k [Myr] 


El 


1 


14 


9.6 


1.9 x 10 6 


1.9 x 10 6 


142 


0.197 


12.57 


E2 





38.2 


10.3 


1.2 x 10 6 


1.2 x 10 6 


129 


0.09 


33.59 


Fl 


4 


9.7 


42 


5x 10 3 


6.5 x 10 5 


2400 


0.047 


2.14 


F2 


9 


28.2 


18.4 


7.3 x 10 5 


7.3 x 10 5 


240 


0.067 


0.54 


F3 


6 


118 


15.4 


2.9 x 10 6 


2.9 x 10 6 


367 


0.1 


0.28 


Gl 


3 


10.14 


45.6 


2.5 x 10 3 


1.1 x 10 6 


4300 


0.011 


0.35 


G2 


3 


13.1 


23.8 


1.5 x 10 4 


6.6 x 10 5 


762 


0.009 


0.009 


G3 


11 


167.4 


92.7 


2.5 x 10 3 


4x 10 6 


7900 


0.1 


44.8 


HI 


4 


11.7 


15.5 


5.6 x 10 5 


1.3 x 10 6 


26 


0.012 


3.65 


H2 


9 


20.1 


17.8 


1.8 x 10 6 


1.8 x 10 6 


360 


0.15 


5.32 


H3 


11 


49.9 


30.2 


1.5 x 10 5 


9.7 x 10 5 


167 


0.28 


9.54 



TABLE 2 

Data for the simulations where the IMBH was captured by a cluster of the CC. The first column shows the ID of the simulation 
(see table[T}. The second column shows the number of stellar collisions triggered by the IMBH. The third column displays the 

TIME OF CAPTURE OF THE IMBH BY A SINGLE CLUSTER. THE FOURTH SHOWS THE DISTANCE FROM THE CENTRE OF THE CC, WHERE THE IMBH WAS 

CAPTURED. The next column gives us the mass of that cluster and the mass of the heaviest cluster in THE CC BY THAT time; i.e. the 

MASS OF THE FORMING UCD. THE SIXTH COLUMN CORRESPONDS TO AN ESTIMATE FOR THE IMBH TO REACH THE CENTRE OF THE CC BY 
DYNAMICAL FRICTION (SEE TEXT). THE LAST TWO COLUMNS SHOW THE TIME THE IMBH HITS A CLUSTER FOR THE FIRST TIME AND THE TIME OF 
THE FIRST STELLAR COLLISION IN THE CC. IN THE PARTICULAR CASE OF SIMULATION Fl THERE WAS A TIDAL DISRUPTION OF A STAR BY THE IMBH. 
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IMBH ID 


Outcome 


Rinit (PC) 


Vinit/Vesc 


T (Myr) 


R|/ n CD (pc) 


M cl (M ) 


# Stellar disr 


1 


capture 


0.0024 


0.89 


14.24 


51.2 


2.5 x 10 3 


1 


2 


capture 


17.32 


0.70 


12.41 


103 


7.85 x 10 5 


2 


3 


satellite 


17.32 


0.71 


34 





9.6 x 10 5 


2 


4 


satellite 


17.32 


0.94 


9.45 


9.5 


5.44 x 10 5 


1 


5 


capture 


17.32 


0.76 


2.35 


118.2 


5 x 10 5 


1 



TABLE 3 

Data for the simulation with five IMBHs in a CC. The first column shows the ID of the IMBH, the second column the outcome of 

THE BH AFTER 35 MYR, WHICH CAN BE EITHER A CAPTURE, A SATELLITE (THE IMBH IS ORBITING A CLUSTER AND WILL EVENTUALLY MERGE WITH 
IT) OR AN ESCAPE (THE IMBH ESCAPES THE WHOLE CC). THE THIRD COLUMN DISPLAYS THE INITIAL DISTANCE OF THE IMBH FROM THE CENTER OF 

THE CC. The fourth corresponds to the initial velocity of the normalised to the local escape velocity from the cluster. The 
fifth gives the time at which the outcome was measured. the sixth shows the final distance of the capturing cluster from the 

most massive cluster of the system, the seed ucd. in this case, imbh #3 is captured by the seed, and thus this distance is zero. 

The seventh is the mass of the capturing cluster at the time of capture. Finally, the last column shows the number of stellar 

collisions in clusters that have been triggered by the imbh. 



IMBH ID 


Outcome 


Rinit (PC) 


Vhiit/Vesc 


T (Myr) 




M cl (M ) 


# Stellar disr 


1 


satellite 


0.0018 


0.56 


1.9 


130.7 


7.5 x 10 5 


2 


2 


satellite 


17.33 


0.97 


23.8 


130.7 


1.6 x 10 6 


2 


3 


satellite 


17.33 


0.99 


8.9 


112.2 


1.44 x 10 6 


3 


4 


satellite 


17.33 


0.89 


9.1 


163.5 


2.2 x 10 6 


1 


5 


escaper 


17.33 


0.56 








5 


6 


escaper 


17.33 


0.59 











7 


escaper 


17.33 


0.88 








1 


8 


satellite 


17.33 


0.90 


10.2 


8.7 


2.7 x 10 6 


1 


9 


satellite 


17.33 


0.56 


5.7 


140.5 


7.75 x 10 6 





10 


satellite 


25.99 


0.72 


4.4 


140.5 


1.47 x 10 6 


1 



TABLE 4 

Same as in table[3]but for ten holes. In this case three IMBHs leave the CC. We find 22 stellar disruptions during the 

SIMULATION. 
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Total mass [xl0 7 M o ] 
1.522 2.9 4.32 5.75 7.18 8.6 11.411 
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Fig. 8. — Outcome of the CC simulations. The x-axis shows the number 
of clusters in each simulation, while the y-axis shows the initial radius of 
the CC. The upper x-axis shows the total mass of the system in Mq . Every 
circle corresponds to a single entry of Table [T] in a way such that the circle 
at the bottom left corresponds to the simulation with ID Al and the circle 
at the top right corresponds to the simulation with ID H4. An open circle 
indicates a simulation where the IMBH finally escaped the CC. On the other 
hand, a filled circle represents a simulation where the IMBH remained bound 
to the system. Next to every circle there are three numbers. The first (black) 
shows the number of clusters hit by the IMBH until either it escapes the CC 
or it is captured by a cluster. The second (red) number is the number of stars 
that are tidally disrupted by the IMBH and the number of star-star collisions 
triggered by the IMBH in the clusters. The third number indicates the initial 
escape velocity at the centre of the CC in kins" 1 . 




80 100 

Time (Myrs) 

Fig. 9. — Cumulative number of IMBH and cluster hits for the simulation 
G3 (inverted, light magenta triangles) and of stellar collisions leading to a 
disruption (blue triangles) as a function of time. 

to reach the centre. For an object with mass m movi ng in 
a system of to tal mass M it is given by (see e.g. Binney & 
|Tremaine|2008] ) 



1.17M r 
InAmVh' 



(5) 



where r is the distance from the centre of the system, V/j is 
the root mean square (RMS) velocity dispersion of the sys- 
tem and In A the Coulomb logarithm, which is of the order of 
unity. From table [2] we can see that in almost half of the cases 
in which the IMBH was retained in the CC, it is captured by 
the most massive cluster of the system, the seed UCD. Tdf is 




10 12 



Time (Myrs) 



Fig. 10. — IMBH speed (dashed, red line with stars) and instantaneous 
escape velocity (solid, green curve with spheres) for the IMBH as a function 
of time in simulation G3. Even though initially the IMBH recoiling speed 
is higher than the required threshold to escape the CC, soon after ~ 0.80 
Myrs the interactions with individual clusters lower its kinetic energy and it 
is trapped in the CC, in the meaning that the speed drops below the threshold. 




, , : 



X(°pc) 

Fig. 11. — Projection in X-Y of all clusters for the simulation in which 
we have initially 5 IMBHs. We show in white, orange, blue, green, yellow, 
and red the clusters that captured the holes (or will capture, if in satellite, see 
text). For clarity we depict all other clusters with the same radius and colour 
(light orange). The green cluster harbours two IMBHs and the blue cluster 
too. The later one merged with an IMBH and after that with another one 
which contained another IMBH. 

in all cases a few tens or hundreds of Myrs. On the other hand, 
when the IMBH gets captured by a smaller cluster (6 out of 
the 1 1 simulations), Tdf is of the order of ~ 1 Gyr, still well 
below a Hubble time. We note also that this analytical calcu- 
lation is an overestimate, because the CC evolves dynamically 
with time and there is a huge accumulation of mass in the in- 
nermost region which will significantly reduce the timescale 
for the IMBH to reach the seed UCD. 

When the IMBH remains bound to the CC, the average time 
for it to hit a cluster is 0.16-0.43Myr. On the other hand, the 
mean time taken by the IMBH to fly through a star cluster is of 
the order of O.lMyr. Hence, after recoiling and before getting 
captured, the IMBH spends 1/3 of its time interacting with 
clusters, so that the possibility to find an IMBH in a cluster of 
a newly formed (less than lOOMyr old) CC is about ~ 30% 

We repeated the exercise with a CC harbouring initially 5 
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Fig. 12. — Same as figure {8} but assuming a fixed loss of kinetic energy 
of 5% after every hit for the fly-throughs. 



Also, reducing the relative kinetic energy loss for fly- 
throughs leads to a reduced number of tidal disruption events, 
but we still find some systems for which the implications are 
similar to the analysis that used a larger loss. 

While the number fraction of IMBH in the mass-range of 
1O 2_4 M in CCs is an unknown, they sink to the centre in a 
time which is much shorter than the Hubble time. The sce- 
nario that we have described here leads to the formation of 
a very massive black hole at the centre of the UCD, with a 
mass that depends on unknowns, such as the formation rate of 
IMBHs in the CC. The internal velocities of th e systems we 
study are not as extreme as those explored by Merritt et al. 
(2009| in the context of hypercompact stellar systems, be- 
cause the seed UCD inherits the central velocity from the re- 
sulting mergers between individual clusters. When the UCD 
is formed, the velocity will roughly be what one can expect 
from a dense stellar system in dynamical equilibrium. A very 
interesting feature of the process of sowing an UCD with an 
IMBH is that independently of whether the IMBH stays in 
the CC or escapes, it triggers star- star collisional disruptions 
in the clusters it hits. This could be envisaged as an electro- 
magnetic signature of the scenario. 



and 10 IMBHs which were distributed with different veloci- 
ties. We find that after some < 30 Myr most of the holes are 
either captured by a single cluster or have formed a hard bi- 
nary with one in regions relatively close to the runaway clus- 
ter, the seed of a UCD which is forming in the CC. We cannot 
follow the further evolution of the system due to the limita- 
tions inherent to our approach. We also note that gas is very 
likely to play an important role in the whole process. In par- 
ticular, in some CCs the oldest cluster is located at the centre 



of the gas cloud (Whit more et al.|2010| ). In our simulations we 
have neglected this, since we are limited by our codes, which 
rely on pure particle dynamics. Still, even if we could actu- 
ally have implemented a (rough) approach for the gas with an 
external force, we believe that the complexity of the problem 
justifies our first approach. We have decided to postpone the 
role of the gas for upcoming work. 
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